Tuning the Mott transition in a Bose-Einstein condensate by multi-photon absorption 
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We study the time-dependent dynamics of a Bose-Einstein condensate trapped in an optical 
lattice. Modeling the system as a Bose-Hubbard model, we show how applying a periodic driving 
field can induce coherent destruction of tunneling. In the low-frequency regime, we obtain the novel 
result that the destruction of tunneling displays extremely sharp peaks when the driving frequency 
is resonant with the depth of the trapping potential ( "multi- photon resonances"), which allows the 
quantum phase transition between the Mott insulator and the superfluid state to be controlled with 
high precision. We further show how the waveform of the field can be chosen to maximize this effect. 
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Recent spectacular progress in trapping cold atomic 
gases [l| has provided a new arena for studying quantum 
many-body physics. In particular, ultracold bosons held 
in optical potentials provide an almost ideal realization 
of the Bose-Hubbard (BH) model :2], in which the model 
parameters can be controlled to high precision. As well 
as their purely theoretical interest, these systems attract 
attention because of their possible application to quan- 
tum information processing |3|. 

The BH model is described by the Hamiltonian 



Hbh = -J^^ \A^ 



H.c. 



U 



^n,(n,~l), (1) 



where ai/a\ are the standard annihilation/creation oper- 
ators for a boson on site i, n^ = a\ai is the number oper- 
ator, J is the tunneling amplitude between neighboring 
sites, and U is the repulsion between a pair of bosons 
occupying the same site. Its physics is governed by the 
competition between the kinetic energy and the Hubbard 
interaction, and thus by the ratio U/ J. When C// J -C 1 
the tunneling dominates, and the ground state of the 
system is a superfluid. As U/ J is increased the system 
passes through a quantum phase transition, and evolves 
into a Mott-insulator (MI) state in which the bosons lo- 
calize on the lattice sites. 

This phase transition was observed experimentally in 
Ref.y by varying the depth of the optical potential. In 
this Letter we propose an alternative method: applying 
an additional oscillatory potential induces coherent de- 
struction of tunneling (CDT), and thus suppresses the 
effect of J. CDT is a quantum interference effect, dis- 
covered in the pioneering work of Ref.|5(, in which the 
period for tunneling between states diverges as their as- 
sociated quasicnergies |g| approach degeneracy. 

Here we show how CDT can be used to control the dy- 
namics of a boson condensate, by means of a novel reso- 
nance effect between U and the frequency of the driving 
field. We consider a one-dimensional BH model, driven 
by a time-periodic potential which varies linearly with 



site number. The Hamiltonian is given by: 



N 



H{t)^HBH+Kf{t)Y, j 



(2) 



where K is the amplitude of the driving field, and f{t) is 
a T-periodic function of unit amplitude that describes its 
waveform. Such time-periodic linear potentials - gener- 
ated by an accelerated lattice for example - have already 
been used in cold-atom experiments 7]. A similar form 
of driving-potential was also recently investigated theo- 
retically 8] in the high-frequency regime (oj > C/), and 
was found to suppress the transition to the superfluid 
regime. 

Here, for the first time, the multi-photon (low- 
frequency) regime is investigated. An unexpected new 
finding is that CDT is now modulated by a set of ex- 
tremely sharp 'resonances' (the contrast between the 
high-frequency behavior and the multi-photon regime in- 
vestigated here is illustrated in Fig^) . This means that 
the Mott-transition can be induced by minute changes in 
experimental parameters. 

Henceforth we put h, = 1 and measure all energies in 
units of J, and set the number of bosons equal to the 
number of lattice sites N. Although the dimension of 
the Hilbert space increases exponentially with iV, in a 
Fock basis H is extremely sparse, with at most {2N — 1) 
non-zero entries per row. Thus despite the rapid increase 
in the dimension of the Hilbert space, this sparsity allows 
us to treat relatively large systems of up to eleven sites, 
and so assess if the effects we observe survive in the ther- 
modynamic limit. 

Our numerical investigation consists of initializing the 
system in the 'ideal' MI state, |^mi) = H'^jlO)' ^^'^ 
evolving the many-particle Schrodinger equation over 
time (typically ten periods of the driving field) using 
a Runge-Kutta method. To study the system's time- 
evolution quantitatively, we measure the overlap of the 
wavefunction with the initial state P{t) — |(^miI*(0)I^- 
For convenience we term the minimum value of P{t) at- 
tained during the time-evolution to be the localization. 
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FIG. 1: (a) The minimum overlap with the Ml-state, or lo- 
calization, reached in a 7-site system with U — 8, during 10 
periods of driving. For oj = 20 (dashed-line) the localization 
peaks at K/lo = 2.4,5.5 - the zeros of Jo(K/uj). When uj 
is reduced to a; = 8 (the first photon resonance, solid-line), 
the peaks become extremely narrow and are centered on the 
zeros of JxiKjio). The diamonds mark the points Kjuj = 3.5 
and 3.8 (see below), (b) Time evolution of the lo = % case for 
three system sizes, 7, 9 and 11 sites. For Kjuji = 3.5, away 
from the resonance, the overlap with the initial state, P(i), 
rapidly drops to zero. For Kjuj — 3.8, at the peak of the res- 
onance, the decay is much slower, indicating that the driving 
field preserves the Ml-state. 



When CDT occurs, the system will remain frozen in the 
Ml-state, and consequently the localization will be close 
to one. Conversely, if the bosons are able to tunnel freely 
from site to site, the value of the localization will be re- 
duced. 

We begin by considering the case of sinusoidal driv- 
ing, f{t) = sinwi. The MI transition is quite soft in 
ID, starting at [/ ~ 4 and developing fully for U > 20. 
Throughout this work we use an intermediate value of 
U = 8. Fig^ shows how the localization in a 7-site 
system varies as the amplitude of the driving field is in- 
creased, while its frequency is held constant at ui — 20. 
For K = the localization has a value of ~ 0.3, demon- 
strating that in the absence of a driving field this value 
of U is indeed insufficient to maintain the Ml-state. Ap- 
plying the driving field causes the localization to steadily 
rise from this value as K is increased from zero, indi- 
cating that the effective tunneling between lattice sites 
is increasingly suppressed, until it peaks at a value close 
to one at K/w = 2.4. As K is increased further, the lo- 
calization goes through a shallow local minimunijbefore 
again peaking at K/to = 5.5. It was observed Q that 
these values of K/u are close to the first two zeros of j7o , 
the zeroth Bessel function. 



Reducing the driving frequency to a lower value, w = 8, 
produces a radically different behavior - the value of the 



localization rapidly drops as K/cu is increased from zero, 
indicating that the field destroys the Ml-state. As K/cu is 
increased further the value of the localization remains ex- 
tremely low except at a series of very sharp peaks. Fig^a 
emphasizes the narrowness of these peaks by showing the 
time-evolution of the system for two values of K. For 
the first, K = 3.5a;, P{t) rapidly falls from its initial 
value, reaching a level near zero within five driving pe- 
riods. There is a small dependence on the system size, 
with the decay occurring more quickly as N is increased. 
At the localization-peak, K = 3.8u;, P{t) decays far more 
slowly with time, so that after twenty periods of driving 
it only falls to a value of ~ 0.9, and only minor depen- 
dence on N is evident. Thus for this value of w, altering 
the amplitude of the field by just 10% produces enormous 
differences in the localization. 

Although the Hamiltonian jSJ is explicitly time- 
dependent, the fact that it is periodic allows us to use 
the Floquet theorem to write solutions of the Schrodinger 
equation as ipit) = exp [— ie^i] 0j(i), where ej is the 
quasienergy, and (j)j (t) is a T-periodic function called the 
Floquet state [g- As the quasienergies are only defined 
up to an arbitrary multiple of w, the quasienergy spec- 
trum possesses a Brillouin zone structure, in precise anal- 
ogy to the quasimomentum in spatially periodic crystals. 
For the Floquet analysis, we work in an extended Hilbert 
space of T-periodic functions Q- In this approach, the 
Floquet states and quasienergies satisfy 



n{t)\<t>,{t))^e,\4>,{t)) , 



(3) 



where 7Y(i) = H{t) — ihd/dt. Working in this extended 
Hilbert space thus reduces the task of calculating the 
time-dependent, driven dynamics of the system to a time- 
independent eigenvalue problem. 

To study the behavior of the quasienergies, we make 
use of a perturbative scheme developed in [lO| to treat 
non-interacting systems, and later generalized in |ll| to 
include interactions. Our procedure is to first find the 
eigensystem of the operator Ho{t) = Ho{t) — ihd/dt, 
where Hq contains terms diagonal in a Fock basis (i.e. 
the driving term and the Hubbard interaction). We are 
then able to use standard Rayleigh-Schrodinger pertur- 
bation theory to evaluate the corrections to this result, 
using the remaining terms of Hbh as the perturbation. 

For the two-site system, a natural basis is given by the 
Fock states {|1,1), |2,0), |0,2)}, where \n,m) denotes 
the state with n bosons on the first site and m on the 
second. Finding the eigensystem of Ho then amounts to 
solving three first-order differential equations, yielding 
the result: 
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-i{U — e+)i -I- i — cosLut 
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0,0, exp 



K 



-i{U — eo)t — i — coswi 



|0_(i)) = (exp[ie_i],O,O). 



(4) 



Imposing the T-periodic boundary condition on these 
states requires setting e_ = and (U — e+/o) = Tnto, 
where m is an integer. Thus in general it is not possible 
to include the full Hubbard- interaction term within Hq, 
depending on its commensurability with uj. To deal with 
this it is necessary to decompose U into a form which du- 
plicates the Brillouin zone structure of the quasienergies 



U 



12 



(5) 



where u is the "reduced interaction", |u| < uj/2. This 
decomposition reveals that only the reduced interaction 
need be included in the perturbation, while the remainder 
of U (an integer multiple of cu) can be retained in Hq. 

To first-order it is easily shown that the three 
quasienergies are given by 



^0 — u and e± = ( u ± 



16^cff 



/2 



(6) 



where the intersite tunneling has been reduced to an ef- 
fective value Jgff = J Jn{K/ui), and n and u are defined 
in Eq[Sl Thus in the high-frequency limit (oj ^ U)^ when 
n — Q and u ^ U , it is clear that the quasienergies cq 
and e+ are degenerate when Jq{K/uj) — 0. In Fig|2Ji we 
show the excellent agreement between the perturbative 
result and the exact quasienergies for a driving field of 
frequency w — 20. Below, in Fig[2l3, we plot the corre- 
sponding value of the localization, and it can be clearly 
seen that the peaks in this quantity are indeed centered 
on the points of closest approach of the quasienergies. It 
may be seen from EqEl that for large values of m, the 
quasienergy separation (e+ — eo) ~ 4 Jgg-/u. The effect of 
u is thus to reduce the amplitude of oscillations in this 
quantity, and so to smear out the avoided crossings of the 
quasienergies. As a result, the peaks in the localization 
are rather broad and overlap each other, and thus the 
localization cannot reach particularly low values. 



Eql^l reveals the particular importance of photon res- 
onances, when U is an integer multiple of the frequency 
of the driving field, U = mo. When this condition is 
satisfied the reduced interaction is zero, and the cross- 
ings between the quasienergies are well-defined. This 
is the origin of the extremely sharp peaks in localiza- 
tion seen in Fig^J^ ior uj — U — %. Away from these 
peaks, the photon-absorption compensates for the en- 
ergy cost of doubly-occupying a lattice site, in analogy to 
the photon-assisted tunneling studied in Rcf. '13|, thereby 
producing low values of localization. In Fig|3; we plot the 
quasienergies for the first photon resonance (ti = 1) for 
the two-site system. For weak fields {K/uj < 2) small 
deviations of the exact quasienergies from the perturba- 
tive result are visible, but for higher field strengths the 
agreement is again excellent. In Fig[2li we plot the local- 
ization produced in this system, and we can note that. 
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FIG. 2: (a) Quasienergy spectrum of a 2-site system, with 
[/ = 8 and ui — 20. Only two of the three quasienergies 
are plotted (the remaining one oscillates weakly about zero), 
which make a series of close approaches to each other as K is 
increased. Red-lines denote the perturbative solutions, which 
agree well with the exact results (black circles), (b) At the 
points of close approach, the tunneling is suppressed and the 
localization peaks. For all field strengths the tunneling is sup- 
pressed with respect to the undriven system, and the localiza- 
tion is thus enhanced, (c) At lower frequencies the behavior of 
the quasienergies changes dramatically. At a; = 8 the system 
is at the first photon resonance {U = ui), and the behavior 
of the quasienergies is described extremely well by the per- 
turbative solutions e — 0,±2J'i{K/u>). (d) As before, the 
localization is peaked at the points of quasienergy crossing, 
but in contrast to (b) the peaks are extremely sharp. 



as seen previously in the 7-site system, the localization 
takes extremely low values except at a set of very nar- 
row peaks. These peaks are precisely aligned with the 
quasienergy crossings at the zeros of J'i{K/uj). 

As we reduce uj still further, we can expect to encounter 
a sequence of higher resonances with similar behavior. 
In Fig|3^ we compare the values of localization produced 
in a 7-site system for the n = 1 and n = 2 resonances. 
The second photon resonance, however, produces a worse 
result than for n = 1. Although sharp peaks are still 
present in the localization, and in agreement with the 
perturbation theory they are indeed centered on the zeros 
of J2{K/lo), the maximum value of localization produced 
is considerably lower. 



This poor localization occurs because not all of the 
quasienergy crossings in the Floquet spectrum of a many- 
site system will occur at precisely the same value of K/uj] 
instead the various crossings occur over an interval. Thus 
although at the peaks in the localization many pairs of 
Floquet states will be degenerate (and tunneling between 
them will be suppressed), other state-pairs will only be 
approximately degenerate and will permit a small, but 
non-zero, degree of tunneling to occur. The major fac- 
tors influencing this effect arise from higher-order terms 
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FIG. 3: Localization produced in a 7-site system (f/ — 8) 
for two forms of periodic driving: (a) sinusoidal, (b) square- 
wave. Both waveforms produce excellent localization at the 
first photon resonance, ui = U, shown by solid-black lines. At 
the second photon resonance (2aj = L'^), shown by red-dashed 
lines, the localization produced by the sinusoidal driving is 
considerably smaller, but the squarewave still produces high 
peaks. 



in the expansion of the single-period time-evolution oper- 
ator [/(T, 0), which manifest as multi-particle tunneling 
and tunneling beyond nearest neighbors. CDT is a quan- 
tum interference effect, which occurs when the dynamical 
phase acquired by a particle in a period of driving pro- 
duces destructive interference, thereby suppressing the 
particle's dynamics. If, however, a 'clump' of ni bosons 
tunnels between sites, the dynamical phase will be rii 
times larger than that for a single boson. Similarly, if 
a boson tunnels between two sites separated by n2 > 1 
lattice spacings, the dynamical phase will be n2 times 
larger. For sinusoidal driving, the single-particle tunnel- 
ing is suppressed when Jn{K/uj) = 0; for these higher- 
order processes also to be suppressed we therefore also 
require J'n{nin2K / lo) ~ for integers rii, 712 = 1, 2 ... A^. 

Clearly this condition cannot be satisfied for sinusoidal 
driving, as the zeros of Jn{x) are not equally spaced. 
Thus to observe good localization properties at high 
photon-resonances we need to to construct a driving field 
f{t) such that the crossings in its Floquet spectrum are 
periodically spaced. This problem was confronted in a 
different context in Ref. 13], where it was shown that such 
a field must be discontinuous at changes of sign. Possibly 
the simplest field of this type, and the most convenient 
for experiment, is a squarewave field. 

In Fig|2b we show the localization obtained in a 7- 
site system driven by a squarewave. Using the same 
perturbative approach as before, it may be shown that 
the quasienergy degeneracies occur for (K/uj) = 2m + 
1 or 2m, depending on whether the order of the res- 
onance, n, is odd or even. Unlike the case of sinu- 
soidal driving, the n = 2 resonance displays good lo- 




FIG. 4: The localization produced in a 5-site system with U = 
8 as a function of the frequency u of the squarewave driving 
field and its amplitude K. When nu = U the localization 
is almost zero (the horizontal blue bands) except at sharply 
defined peaks (red). Between the bands localization is good 
(red). The n = 1,2,3 and 4 resonances are marked on the 
right. The upper-left triangle (blue) displays poor localization 
and little structure, corresponding to the non-perturbative 
regime. 



calization, comparable to that obtained for n = I. A 
contour plot showing the localization as a function of 
both K/lu and w"""^ is presented in Fig^] The prominent 
horizontal bands correspond to the photon resonances 
{uj^^ — n/U), which are punctuated by a series of nar- 
row peaks at which the localization is preserved. This 
plot also clearly shows the division between the fairly fea- 
tureless, poorly-localized, "weak-driving" regime to the 
upper left, and the "strong-driving" regime which shows 
the resonance features. For the latter, the dynamics of 
the system are dominated by the combined effect of the 
driving field and the Hubbard interaction, and thus is 
well-described by our form of perturbation theory. Fig0] 
allows us to locate the boundary between the two regimes 
quite accurately as K/lu ~ {2U)/lo. 



In summary, we have investigated the dynamics of the 
BH model under a periodic driving field. For high fre- 
quencies I3] the field can be used to inhibit tunneling by 
means of CDT and thus stabilize the Ml-state. Lower- 
ing the frequency, however, reveals the existence of res- 
onance effects which can be used to selectively destroy 
or preserve the Ml-state. Lower driving frequencies have 
the added advantages in experiment that they heat the 
condensate less, and will not drive transitions to higher 
Bloch bands thereby invalidating the single band model. 
The extremely narrow width of the resonance features 
indicates that it should be possible to control the Mott 
transition very precisely in this manner. 
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